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ABSTRACT 

We use 62,185 quasars from the Sloan Digital Sky Survey DR5 sample to explore 
the relationship between black hole mass and luminosity. Black hole masses were 
estimated based on the widths of their H/3, Mgll, and ClV lines and adjacent continuum 
luminosities using standard virial mass estimate scaling laws. We find that, over the 
range 0.2 < z < 4.0, the most luminous low-mass quasars are at their Eddington 
luminosity, but the most luminous high-mass quasars in each redshift bin fall short of 
their Eddington luminosities, with the shortfall of order ten or more at 0.2 < z < 0.6. 
We examine several potential sources of measurement uncertainty or bias and show 
that none of them can account for this effect. We also show the statistical uncertainty 
in virial mass estimation to have an upper bound of ~ 0.15 dex, smaller than the 0.4 
dex previously reported. We also examine the highest-mass quasars in every redshift 
bin in an effort to learn more about quasars that are about to cease their luminous 
accretion. We conclude that the quasar mass-luminosity locus contains a number of 
new puzzles that must be explained theoretically. 

Key words: black hole physics — galaxies: evolution — galaxies: nuclei — quasars: 
general — accretion, accretion discs 



1 INTRODUCTION 

Supermassive black holes (SMBH), with masses between 
~ 1O 6 M0 and ~ 1O 9 M0, are found at the center of nearly 
every galaxy where there have been sensitive searches. While 
we suspect that the seeds for these SMBH might all have a 
common origin, their formation mechanism is not well un- 
derstood. Many galaxies at redshifts z ~ 2 contain quasars, 
i.e., SMBH in the midst of luminous accretion. The Soltan 
argument (Merritt & Ferrarese 2001 ) suggests black hole 
masses are largely accounted for via growth due to lumi- 
nous accretion. There are far fewer quasars at low red- 
shift l|Schmidt fc Green 1983llRichards et al. 2006a)) . imply- 
ing that at some point, SMBH cease their luminous ac- 
cretion. The quasar turnoff mechanism is not well under- 
stood (Thacker et al. 2006). Finally, the black hole mass - 
stellar velocity (M — a) relation ( Ferrarese & Merritt 2000 
IGebha rdt et al. 2000 ) suggests that SMBH are in some way 
co-evolving with their host galaxies, but the M — a rela- 
tion merely describes an end state and is not a theoretical 
explanation. 

As we know of no rapid process by which SMBH can lose 
mass, the evolutionary tracks for SMBH consist of stages 
at progressively higher masses. These stages must, in or- 
der, involve (1) a formation mechanism, (2) a period of lu- 
minous growth (the 'quasar phase'), perhaps along with a 



period of nonluminous growth ('turnoff'), and (3) a period 
in which SMBH lie at the centers of galaxies without sub- 
stantial growth, as we observe them today. Much current 
theoretical research concerns the origin and turnoff phases 
of quasar evolution, while the 'quasar phase' appears to be 
relatively well understood. 

In this paper we find a new feature of SMBH evolu- 
tion during their quasar phase. We examine the evolution 
of the quasar locus in mass-luminosity space as a function 
of redshift. The M — L locus is traditionally shown with all 
quasars on the same plot (as in Figure [I). The large size of 
the Sloan Digital Sky Survey (SDSS) DR5 quasar catalogue 
l|Schneider et al. 2007P allows a subdivision into several red- 
shift bins, each containing thousands of quasars. 

The Eddington limit produces an absolute upper bound 
on the luminosity of quasars which is proportional to the 
black hole mass, L E dd = 1.3 x 10 46 (M/10 8 M Q ) erg s^ 1 
(Shapiro & Tcukols ky 1983[ ). While strictly applicable only 
for a spherical accretion flow of ionized gas, models for more 
realistic accretion configurations with rotation are still lim- 
ited by a luminosity of this order (except under special cir- 
cumstances, e.g. Begelman 2002). Using a smaller data set 
(N = 733) than SDSS, Kollmeier et al. (2006) appeared to 
confirm the applicability of L E dd to quasars, using an M — L 
locus (similar to Figure [1} to show that the most luminous 
quasars reach but do not exceed L E d- For the lowest black- 
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hole masses at every redshift, we confirm the conclusion of 
l|Kollmeier et al. 2006)1 . but we show that the quasars with 
the highest SMBH mass at every redshift fall well short of 

LEdd- 

In § [2] we review the methods used to estimate masses 
and bolometric luminosities. While we have added nothing 
original to this methodology, the remainder of our results 
are entirely dependent upon its accuracy. In § [3] we subdi- 
vide the SDSS DR5 quasar catalogue by redshift and show 
that the quasar mass-luminosity distribution does not match 
what we should expect given our current theoretical under- 
standing of quasar accretion. In particular, we show that 
instead a sub-Eddington boundary (SEB) is present in each 
redshift bin. We also use our limited statistics to make a first 
estimate for how the SEB evolves with redshift. In § [H we 
consider several alternative explanations for the discrepan- 
cies between the observed quasar locus and the Eddington 
luminosity, focusing on potential sources of measurement 
uncertainty or bias. In § [5] we comment on the potential 
implications of our new results. 



2 VIRIAL MASS ESTIMATION 

The remainder of this work relies upon an ability to ac- 
curately estimate the bolometric luminosities and central 
black hole masses of quasars at cosmological redshifts. 
The primary results in this work are drawn from the 
Shen et al. (2008) virial mass catalogue for SDSS DR5 
IjSchne idcr et al. 2007) quasars. The luminosity determina- 
tion is fairly straightforward and uses the relatively settled 
techniques discussed in Richards et al. (2006). The mass 
determination, on the other hand, uses relatively new tech- 
niques, and has a larger uncertainty. Therefore, we review 
here the basic assumptions in virial mass estimation. Poten- 
tial sources of error or bias in the Shen et al. virial masses 
are discussed in more detail in § [4] 

The determination of black hole masses from spectral 
emission lines relies upon two basic assumptions: (1) that 
the orbital velocity of gas in the broad-line region (BLR) 
is dominated by the virial velocity due to the central black 
hole and (2) that there is a scaling relationship between 
luminosity and radius. As a result of (1), we can calculate 
the mass Mbh of the black hole from emission lines in the 
BLR as 



Table 1 . Summary of the virial mass estimates used 



Mbi 



^BLR^BLR 
G 



(1) 



where A^gLR, ' s ^ ne ra dius and vgLR, the velocity of gas 
emitting the BLR spectral lines. Marconi et al. (2009) sug- 
gest corrections to the virial approximation for radiation 
pressure might be needed at high luminosity, particularly 
when using CIV lines to determine vg^F!/ This first scal- 
ing relationship is based upon black hole masses determined 
using reverberation mapping, which uses the time delay be- 
tween variability in the continuum and emission lines to de- 
termine ^blj^ (cf. Peterson & Home 2004). 

Virial mass estimates also require assumption (2), an 
empirical scaling relationship very close to the LaR 2 that we 
would expect for a black body l|Bentz et al. 2009")) . although 
thermal processes in the accretion disc are likely to be sub- 
stantially more complex. This second assumption transforms 
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Redshift 


Source 


H/3 mass 


6.91 


5100 


0.50 


0-0.872 


VP06 


Mgll mass 


6.51 


3000 


0.50 


0.393 - 2.252 


MD04 


CIV mass 


6.66 


1350 


0.53 


1.518-4.875 


VP06 



|T]) into a scaling relation of the form 
log(M/M©) = A 



(2) 



+ log 



( FWHM(H/3) \ 2 / XLx(B A) x ' 
\^ 1000 km/s J \ 10 44 erg/s 



where FWHM is the Full Width at Half-Maximum of the 
corresponding line profile, L\ is the luminosity per unit 
wavelength at rest-frame wavelength A = B, and A, B, and 
C are constants. Vestergaard & Peterson (2006) determined 
these constants in black hole mass scaling relations based on 
H/3 and CIV emission lines (VP06), while McLure & Jarvis 
(2002) and McLure & Dunlop (2004) developed a mass re- 
lation (MD04) by scaling Mgll-based estimates against H/3- 
based estimates. The scaling relations used in this work are 
summarized in Table [1] Multiple scaling relations are re- 
quired because SDSS spectra only cover the range 3900-9100 
A, so that none of these three BLR emission lines are acces- 
sible over the entire redshift range of the SDSS catalogue. 

The additive constants A are determined by calibrat- 
ing virial mass scaling relationships against other meth- 
ods for black hole mass estimation as part of a 'black 
hole mass ladder' ( Pet erson fc Home 2004). Direct estimates 
based upon local stellar and gas kinematics (cf. Ferrarese 
& Ford 2005 provide our best estimates for nearby black 
holes. Reverberation mapping is calibrated against these es- 
timates, and the virial mass scaling relations are in turn 
calibrated against reverberation mapping. These calibra- 
tions use very few quasars compared to the 62,185 in 
the Shen et al. (2008) sample: only 28 quasars are used 
to calibrate the H/3 and 27 to calibrate the CIV scaling 
relations (Vestergaard & Peterson 2006). 

There is a small overlap in redshift between the H/3 
mass and Mgll mass samples as well as one between the 
Mgll mass and CIV mass samples (Table [T}. Shen et al. 
(2008) show that the agreement between mass estimates for 
the same object taken using the H/3 and Mgll relations (0.22 
dex dispersion), is better than between estimates using the 
Mgll and CIV relations (0.34 dex, and correlated with the 
ClV-Mgll blueshift). The CIV mass estimates may be less 
accurate and we therefore focus here almost exclusively on 
masses obtained using the other two scaling relations. We 
consider the implications of this disagreement further in § [4] 



3 THE MASS-LUMINOSITY RELATION 

Black hole masses for 62,185 of the 77,429 SDSS DR5 
quasars were determined by Shen et al. (2008) using the scal- 
ing relations summarized in Table [T] 10,605 in the H/3-mass 
sample, 42,035 in the Mgll-mass sample, and 14,565 in the 
ClV-mass sample. The catalogue includes 3505 quasars with 
both H/3 and Mgll masses and 3427 objects with both Mgll 
and CIV masses. Because the ClV-based mass estimates may 
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Figure 1. The quasar locus in the mass-luminosity plane for all 
quasars from 0.2 < z < 2.0, using virial masses estimated by Shen 
et al. (2008) with H/3 and Mgll lines and bolometric luminosities 
using the techniques of Richards et al. (2006) The dashed line is 
drawn at the Eddington luminosity as a function of mass. The 
colour indicates whether the quasar is at 0.2 < z < 0.8 (red), 
0.8 < z < 1.4 (yellow), or 1.4 < z < 2.0 (green). 

be less accurate, our main results are derived from 49,135 
DR5 quasars in the H/3 and Mgll mass samples. 

Figure [1] displays virial mass estimates and bolomet- 
ric luminosities for all quasars in the H/3 and Mgll mass 
samples. The most striking feature is that the quasar locus 
seems bounded by Lsdd (dashed line), as already shown by 
Kollmeier et al. (2006). However, while the Lsdd bound is 
tight at most masses, we also note a slight departure from 
the Lsdd bound for M > 10 9 Mq. It is also apparent that 
the locus of quasars at 1.4 < z < 2.0 (green) is different 
than the locus at 0.2 < z < 0.8 (red). A proper investi- 
gation of this possible 'sub-Eddington boundary' (hereafter 
SEB) therefore requires a subdivision of Figure [T] into red- 
shift bins. 



3.1 Quasar mass evolution 

We have therefore divided the H/3 mass and Mgll mass sam- 
ples into 10 redshift bins of size 0.2 in z. Table [2] contains 
summary statistics on the objects in each bin. Within each 
bin, there is a distribution of black hole masses spanning 
~ 2 dex, as shown in Figure [2] 

It has been known for a long time that quasars are 
'downsizing', or that the brightest quasars at higher-redshift 
are more intrinsically luminous than the brightest quasars at 
lower-redshift (|Schmidt 1968|) . Similarly, quasars above the 
peak in each of the higher-redshift mass distributions in Fig- 
ure[3]he at masses with substantially smaller populations in 
the quasar mass distributions at lower redshift. This shows 
that many higher-mass quasars turn off by lower redshift 
(and disappear from the sample) rather than become less 
luminous (but remain in the sample at lower luminosity). 
We discuss quasar turnoff in more detail in Paper II. While 



Table 2. Summary statistics on quasars in the 10 redshift and 
emission line bins 



ID z N <logL> a L < log M/M > a M 
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0.2-0.4 


2690 


45.25 


0.20 


8.27 


0.44 


2 


0.4-0.6 


4250 


45.54 


0.25 


8.44 


0.42 


3 


0.6-0.8 


3665 


45.89 


0.25 


8.69 


0.39 


Mgll 














4 


0.6-0.8 


4727 


45.80 


0.29 


8.59 


0.32 


5 


0.8-1.0 


5197 


46.02 


0.30 


8.76 


0.31 


6 


1.0-1.2 


6054 


46.21 


0.26 


8.89 


0.29 


7 


1.2-1.4 


7005 


46.32 


0.27 


8.96 


0.29 


8 


1.4-1.6 


7513 


46.43 


0.27 


9.07 


0.28 


9 


1.6-1.8 


6639 


46.57 


0.24 


9.18 


0.29 


10 


1.8-2.0 


4900 


46.71 


0.22 


9.29 


0.30 




8 9 



Log (M/M ) 

Figure 2. The mass distribution for three different redshift 
bins within our sample. The red (0.2 < z < 0.4) and green 
(0.6 < z < 0.8) lines are from the H/3 mass sample and the purple 
(1.2 < z < 1.4) line is from the Mgll mass sample. These mass 
distributions are not completeness-corrected, and the low-mass 
tails of these distributions are likely lowered by SDSS magnitude 
selection. However, high-mass quasars can be detected even at 
low L^ddi so a lack of high- mass quasars at low redshifts cannot 
be ascribed to selection. 



the low-mass quasar distribution also varies with redshift, 
the SDSS detection limit is at a fixed magnitude. Therefore, 
lower-mass and more distant quasars must be closer to their 
Eddington luminosity in order to be bright enough to be in- 
clude in the SDSS catalogue. As a result, the low-mass tails 
of these mass distributions may be skewed by SDSS selec- 
tion. However, high-mass quasars can be detected even at 
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Figure 3. The SDSS quasar locus of the H/3 mass sample in 
the M — L plane at redshift 0.2 < z < 0.4. The locus should be 
bounded by SDSS detection limits, L^dd (dashed line), and 
on the high-mass end by an unknown mechanism responsible 
for quasar turnoff. In practice, there appears to be an addi- 
tional sub-Eddington boundary with slope below that of LEdd- 
The bright-object SDSS saturation limit does not intersect the 
quasar locus. 



low LEdd, and a lack of high- mass quasars at low redshift 
cannot be ascribed to selection. 
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Figure 4. The quasar number density for the quasars shown in 
Figure [3] at 0.2 < z < 0.4 as a function of luminosity in seven 
different mass slices, each of width 0.25 dex in log M/Mq: black 
(7.25-7.5), red (7.5-7.75), yellow (7.75-8.0), green (8.0-8.25), cyan 
(8.25-8.5), blue (8.5-8.75), and purple (8.75-9.0). The dashed line 
is drawn proportional to L~ 2 and is normalized to the purple 
(8.75-9.0) curve. 



3.2 The Mass-Luminosity Plane at 0.2 < z < 0.4 

The brightest quasars are more intrinsically luminous at 
higher redshift, and similarly Figure [2] demonstrates that 
the biggest central black holes are more massive at higher 
redshift. Further, Figure [1] demonstrates that the most lu- 
minous quasars at every mass are near LEdd- It is there- 
fore natural to believe that quasar luminosity downsizing 
and quasar mass downsizing are simultaneous, such that the 
most massive and most luminous quasars decline with equal 
speed in both mass and luminosity towards lower redshift, 
remaining near Lsdd at every redshift. 

In Figure[31 we present the quasar locus at 0.2 < z < 0.4 
in the M — L plane. In some papers this plane is plotted 
with luminosity on the abscissa. We prefer to put mass on 
the abscissa as mass is a less variable property of the ob- 
ject. The origin of the boundaries of this locus are mostly 
understood: The SDSS DR5 selection has magnitude lim- 
its due to detector sensitivity at i ~ 22 (which we have 
labeled as Detection limit in Figure [3]) and SDSS satu- 
ration at i < 16 (which does not bound any of our quasar 
loci). There is also a high-mass limit. As larger SMBH do 
exist at higher redshift, we have labeled this limit as the 
Quasar turnoff. This limit is discussed in detail in Paper 
II (|Steinhardt fc Elvis 2009p . The dashed line is drawn at 
LEdd(M). As in Figure [TJ there are no quasars statistically 
exceeding L E dd- 

More strikingly, there is a tighter bound than Lsdd on 
the maximum luminosity at masses M > 10 7 ' 5 M Q . We term 
this the Sub-Eddington Boundary (SEB; the red line in Fig- 
ure [3]). This SEB is much more prominent than in the entire 
sample (Figure^, implying a redshift evolution of the SEB. 
To determine the shape of the SEB, we consider the quasar 



Table 3. Best-fitting exponential decays N on L k for the quasar 
luminosity function in different mass bins at 0.2 < z < 0.4 



log M/Mq 


slope k 


x 2 /dof 


7.25-7.5 


3.69 ±0.33 


0.41 


7.5-7.75 


4.69 ±0.59 


3.23 


7.75-8.0 


2.96 ±0.35 


1.80 


8.0-8.25 


2.01 ±0.22 


4.77 


8.25-8.5 


2.25 ±0.22 


3.74 


8.5-8.75 


1.79 ±0.15 


3.62 


8.75-9.0 


1.96 ±0.17 


1.01 



luminosity distribution as a function of mass. Figure|4]shows 
the quasar number density as a function of luminosity in dif- 
ferent mass bins including Poisson errors. The distributions 
in each mass bin have a peak number density with a de- 
cline at high luminosity (see § [4} described in Table [3] At 
high mass, where a wide range of luminosity is visible above 
the peak number density and the maximum luminosity is 
sub-Eddington, the best-fitting exponential decay is ~ L~ 2 . 
At low mass, where SDSS detection limits the sample to 
a narrow range of luminosity and the most luminous ob- 
jects lie near LEdd, the quasar number density decline may 
be steeper with increasing luminosity than at high mass. A 
majority of the best-fitting declines have x 2 /DOF > 3, sug- 
gesting that the falloff may not be purely exponential. Since 
the quasars at lowest mass only allow a fit from 4-5 points, 
it is also possible that these declines are all ~ L~ 2 for much 
of their luminosity range but steeper near the SEB. 

The peaks in each mass bin do not lie at the SDSS low- 
magnitude cutoff, but rather are true peaks in the luminosity 
number density. The quasar luminosity function declines as 
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Figure 5. The SEB as approximated by the 95th percentile lu- 
minosity above peak number density at 0.2 < z < 0.4. The best- 
fitting linear approximation has slope is a = 0.37 ± 0.02, well 
below the a = 1 slope of L^dd (white, alternating dashes). The 
indicated uncertainties are derived via bootstrapping and the fit 
has a very low x 2 /DOF of 0.06. Note that the SEB at every mass 
lies well away from the SDSS detection and saturation limits (long 
and short dashed lines, respectively). 



Figure 6. The quasar distribution in the mass-luminosity plane 
at 0.6 < z < 0.8 as measured using H/3 masses. The sub- 
Eddington boundary (red, dashed) is fit from 95th percentile 
points (red) above the peak number density (blue, dashed) and is 
different than that at 0.2 < z < 0.4 (Figure(3)l but is still present. 
The black dashed line is drawn at a bolometric luminosity ap- 
proximately corresponding to i = 19.2 for a typical quasar SED 
at this redshift. 



a power law in luminosity ~ L~ 2 (cf. Richards et al. 2006b, 
Amarie et al. 2009) at luminosities above the peak number 
density. For comparison, a decline proportional to L~ 2 is 
shown in Figure 2] Since the ~ 0.4 dex mass uncertainty 
is large compared to the 0.25 dex mass bins, a majority 
of quasars should randomly lie in the wrong mass bin and 
~ 1/3 of quasars would not even lie in a bin adjacent to their 
correct mass bin. Detailed fit parameters for the decline are 
likely not credible and have a high x 2 /DOF as in Table [31 
but a power-law decline ~ L~ 2 is consistent with the data 
in each of the higher mass bins where quasars take on the 
widest range of luminosity. 

We define the luminosity boundary L au toff(M) ro- 
bustly in each mass bin as the 95th percentile quasar lu- 
minosity above the peak. In each mass bin, the uncertainty 
can be estimated using bootstrapping (repeatedly choosing 
a set of N quasars randomly with replacement from the N 
quasars in each mass bin) , finding the standard deviation of 
the resulting 95th percentile luminosities. Figure [5j shows 
Lcutoff(M) for the 0.2 < z < 0.4 redshift bin. The best 
straight-line fit for L cuto ff has slope a — 0.37 ± 0.02, well 
below OLEdd = 1 (Figure [3]l . An SEB is strongly required. 
The low x 2 /DOF of 0.06 for this linear fit is likely evidence 
of correlated uncertainties, perhaps due to the large mass 
uncertainty placing individual objects in incorrect bins. As 
shown in § 13.31 x 2 /DOF for this redshift bin is atypically 
low. 

The highest-mass bins lie further from our best-fitting 
lines but also have larger uncertainties due to a lower quasar 
number density, while quasars with extremal mass estimates 
are most likely to have been placed in the wrong bin. It is 
possible that the SEB would be well-fit by two linear com- 



ponents given lower-uncertainty measurements of M and L. 
We consider whether such measurement uncertainties might 
be responsible for the a < 1 slope of the SEB in § [4] 

3.3 Evolution of the sub-Eddington boundary 

Figure [H] shows quasar loci in the L — M plane for each 
of the 12 redshift bins as contour plots. An SEB is de- 
tected in each panel, although the location and slope ap- 
pear to evolve with redshift. An offset at which quasars 
fall short of Lsdd when using Mgll masses has been 
previously reported (Kollmci er et al. 20061 [Shen et al. 20081 
|Gavignaud et al. 2008||Trump et al. 20090 . The effect found 
here is different as the offset is mass dependent and changes 
slope with redshift (Figures [6] : the maximum Eddington 
ratio at higher masses is below that at lower masses. 

While the three CIV mass samples have known flaws, 
they are included in Figure [8] in order to demonstrate that 
the best available evidence suggests that the SEB continues 
to exist at redshifts higher than the z ~ 2.0 limitations of 
SDSS spectra-based Mgll mass estimation. 

In § 13.21 we showed that the 0.2 < z < 0.4 SEB appears 
to be well- fit by a linear L = aM + Lq, as in Figure [3] We 
therefore parametrize the SEB in the same manner in each 
redshift bin. In Table [4] we show the results of these fits for 
0.2 < z < 2.0. At every redshift, the SEB fit parameters 
shown in Table [3] show a slope at least 3.7a below that of 
the Eddington luminosity. 

These best-fitting parameters depend upon the peak 
number density. At most combinations of mass and redshift, 
the peak number density lies at a luminosity higher than 
many objects in the SDSS catalog. However, while the SDSS 
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Figure 8. A contour plot of the mass-luminosity distribution in 12 different redshift ranges with the best-fitting SEB (red, dashed), 
fn each redshift bin, the quasar number density has been normalized to the peak number density. The top three panels are from our 
low-redshift sample with H/3-based mass estimates, the middle six panels are from our medium-redshift sample with Mglf-based mass 
estimates, and the bottom three panels are from our high-redshift sample and use CIV-based mass estimates. The three redshift intervals 
with CIV-based masses suffer from suspected flaws in mass estimation (Shen et al. 2008); however, they are included to suggest, with 
the best available data, that the SEB continues to higher rcdshifts. 



catalog is nearly complete for objects brigher than i = 19.2, 
fainter objects are only included if they are "serendipitously" 
selected as a ROSAT source, FIRST source, etc. In Paper II, 
we compare this serendipitous sample more closely with the 
remainder of the SDSS catalog and consider their views of 



the low- luminosity end of the mass- luminosity quasar distri- 
bution. While i = 19.2 does not translate directly to a bolo- 
metric luminosity because the quasar spectal energy distri- 
bution can vary, a bolometric luminosity typical of quasars 
in the SDSS catalog with i — 19.2 is indicated in Figures [5] 
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Figure 7. The quasar distribution in the mass-luminosity plane 
at 1.0 < z < 1.2 as measured using Mgll masses. The sub- 
Eddington boundary (red, dashed) is fit from 95th percentile 
points (red) above the peak number density (blue, dashed) and 
has a slope much closer to the Eddington limit than it does at 
0.2 < z < 0.4 (Figure l3t but is still not parallel. The quasar pop- 
ulation approaches L^dd a t lower mass but does not at higher 
mass. The black dashed line is drawn at a bolometric luminosity 
approximately corresponding to i = 19.2 for a typical quasar SED 
at this redshift. 



Table 4. Fit parameters at different redshift for the maximum 
luminosity L = aM + Lq. This is a two-parameter fit to typically 
8-10 total points. Both the Hp and Mgll views at 0.6 < z < 0.8 



are included. 


Redshift 


a 




Lo 




x 2 /dof 


H/9 












0.2-0.4 


0.37 


0.02 


42.63 


0.18 


0.06 


0.4-0.6 


0.45 


0.03 


42.13 


0.22 


0.28 


0.6-0.8 


0.60 


0.06 


41.07 


0.49 


0.56 


Mgll 












0.6-0.8 


0.61 


0.10 


41.03 


0.83 


1.11 


0.8-1.0 


0.67 


0.09 


40.60 


0.79 


0.90 


1.0-1.2 


0.67 


0.05 


40.71 


0.48 


0.85 


1.2-1.4 


0.73 


0.05 


40.24 


0.48 


1.13 


1.4-1.6 


0.68 


0.08 


40.66 


0.72 


1.10 


1.6-1.8 


0.50 


0.10 


42.35 


0.89 


1.64 


1.8-2.0 


0.42 


0.06 


43.20 


0.51 


1.06 



and [7] It is possible that peak luminosity has been overesti- 
mated because of incompleteness below i = 19.2. The entire 
luminosity distribution at fixed mass and redshift typically 
spans between 0.8 and 1.5 dex (the implications of this are 
discussed further in Paper II). In Tabled] we consider the 
possibility that the peaks are overestimated due to selection 
and show the best-fitting SEB as defined using a peak 0.2 
dex lower in luminosity at every mass and redshift. Because 
of the sharp decline in number density near the SEB, an 0.2 
dex shift in the peak corresponds to a smaller shift in the 
95th percentile object used to estimate the SEB. We also 
consider shifting only peaks at the lowest 1.0 dex of mass at 



Table 5. Fit parameters at different redshift for the maximum 
luminosity L = aM + Lo using the 95th percentile objects above 
or within 0.2 dex of the peak number density at each mass and 
redshift. Best-fitting parameters using the SDSS peak at high 
mass but this shifted peak at low mass are also considered in an 
attempt to bias the SEB determination as far as possible towards 
a slope of 1. 



Redshift 


Original Slope 


Slope, 0.2 below Peak 


Slope, biased peaks 


H/3 








0.2-0.4 


0.37 ±0.02 


0.34 ±0.03 


0.37 ±0.02 


0.4-0.6 


0.45 ± 0.03 


0.48 ± 0.02 


0.49 ±0.02 


0.6-0.8 


0.60 ±0.06 


0.61 ± 0.06 


0.64 ± 0.06 


Mgll 








0.6-0.8 


0.61 ±0.10 


0.60 ±0.08 


0.64 ±0.10 


0.8-1.0 


0.67 ±0.09 


0.65 ± 0.06 


0.69 ±0.07 


1.0-1.2 


0.67 ±0.05 


0.66 ± 0.06 


0.68 ±0.06 


1.2-1.4 


0.73 ±0.05 


0.72 ± 0.05 


0.74 ±0.06 


1.4-1.6 


0.68 ±0.08 


0.67 ±0.08 


0.70 ±0.08 


1.6-1.8 


0.50 ±0.10 


0.51 ± 0.09 


0.54 ±0.10 


1.8-2.0 


0.42 ± 0.06 


0.42 ± 0.06 


0.45 ±0.06 




Redshift 

Figure 9. Evolution of the best-fitting SEB slope displayed in 
Table [4] The green measurements are from the H/3 mass sample, 
while blue points are from the Mgll mass sample. The best-fitting 
linear evolution for the SEB slope is a poor fit, with x 2 /DOF of 
7.84. The magenta points use possible corrections to Mgll masses 
discussed in § 14.51 



each redshift in an attempt to bias the SEB determination 
as far as possible towards a slope of 1. 

In Figure [9] we show the redshift evolution of the SEB 
slope a,. All of the a lie within a narrow range of standard 
deviation 0.12. The best linear fit to the slope as a function of 
redshift is a(z) = (0.12 ± 0.08)z + (0.41 ± 0.09), as in Figure 
El but is a poor fit, with x 2 /DOF of 7.84. This analysis 
cannot exclude the possibility that the SEB takes on a slope 
independent of redshift. 
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Log L /L Ed d 

Figure 10. The distribution of Eddington ratios in five mass 
bins of width 0.25 dex in log M/Mq at 0.2 < z < 0.4 from Figure 
H black (7.75-8.0), red (8.0-8.25), cyan (8.25-8.5), blue (8.5-8.75), 
and purple (8.75-9.0). The low-Eddington ratio boundary is likely 
due to SDSS magnitude limitations. We use a KS test to deter- 
mine whether the overlapping portions of these distributions that 
would pass SDSS selection are identical. Distributions have been 
normalized to 1.0 at peak, and only the portions of each distribu- 
tion where the apparent magnitude corresponding to L is above 
the SDSS detection threshold are shown. 

4 IS THE SUB-EDDINGTON BOUNDARY 
DUE TO MEASUREMENT ERROR? 

We can divide the explanations for the SEB into four possi- 
bilities: 

(i) SDSS selection excluding high-luminosity, high-mass 
quasars from the DR3 and DR5 catalogues. 

(ii) Measurement errors resulting in either an underesti- 
mated bolometric luminosity or incorrect fit parameters for 
the spectral lines used to estimate M. 

(iii) Incorrect virial mass scaling relations for higher-mass 
quasars. 

(iv) Physical effects limiting luminous quasar accretion 
more strongly than the Eddington limit. 

We consider the first three possibilities in this section, with 
a summary of the potential explanations considered in Table 
[(J We briefly consider the implications of these boundaries 
being physical in §0 A full consideration of possible physical 
causes for the SEB is beyond the scope of this paper. 

4.1 Statistical significance 

We first consider whether the SEB is statistically significant, 
or whether the L/LEdd distribution is consistent with being 
identical in different mass bins. Figure [TD] presents the dis- 
tribution of quasar Eddington ratios at 0.2 < z < 0.4 in the 
five most populous mass bins (10 775 < M/Mq < 10 9 ). 
The SDSS detection limit prevents lower-mass bins from 
containing objects at lower Eddington ratios than shown. 
The distributions show a clear trend to higher L/LEdd at 




Log L/L Edd (arbitrary units) 



Figure 11. The five distributions from Figure llOl translated to 
match their median Eddington luminosites. Distributions have 
been normalized to their peak. 



lower masses. This trend is sufficiently strong that there is 
almost no overlap between bins of very disparate mass. That 
these L/LEdd distributions are distinct can be shown via a 
comparison the pair of bins at AlogM greater than the 0.4 
dex mass uncertainty with the largest number of quasars in 
the overlapping region. We use a KS test to quantify whether 
the overlapping portions of a pair of distributions are consis- 
tent with being drawn from the same distribution. For the 
red (8.0-8.25 in log M/Mq) and blue (8.5-8.75) distributions 
from Figure [TU] the KS test yields a D value of 0.5862, and 
the probability that these are drawn from the same distri- 
bution is 3.6 x 10 -50 , i.e., with > 49a confidence these are 
statistically different distributions. Thus, the apparent SEB 
are not merely artifacts of small number statistics. The dif- 
ference in the medians of this pair of distributions separated 
by 0.5 dex in M is 0.34 dex in L/LEdd, or 0.16 dex in L. The 
quasar luminosity distribution at 0.2 < z < 0.4 is neither 
mass-independent nor linear in mass, but rather has some 
sub-linear dependence. 

The luminosity distributions in each mass bin appear 
to have similar shapes (Figure [4}. To test this, the L/LEdd 
distributions are shifted to a common median value. The 
KS test on the shifted red (8.0-8.25) and blue (8.5-8.75) 
distributions now yields a D value of just 0.0747, and the 
corresponding probability that these two distributions have 
the same shape is 0.358. Translating all five mass bins to 
a common median (Figure QT} confirms that these L/LEdd 
distributions are all similar. This similarity is particularly 
striking because the total (sum at all masses) quasar lu- 
minosity function has a similar shape at different redshift 
URichards et al. 2006a; Amari e~et al. 2009]l . These similari- 
ties might hint at an underlying cosmic structure for quasar 
populations spanning a wide range of mass and redshift. 

We conclude that the SEB indeed sub-Eddington with 
very high statistical certainty. Similarities in the shape of 
luminosity distributions at different black hole masses might 
provide a useful hint as to the origin of the SEB. 
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4.2 SDSS selection 

SDSS does not generate spectroscopic data on every object 
in the survey, but rather only on those objects (including all 
quasar candidates) designated as worthy of followup based 
upon their photometry and observed apparent magnitudes 
in five spectral bands ugriz ((Richards ct al. 2002 ) . Every ob- 
ject identified as a quasar candidate is entered into the spec- 
troscopic queue. The DR5 SDSS spectroscopic footprint is a 
subset of the DR5 photometric footprint. Regardless of the 
original photometric classification, any object identified as a 
quasar from its spectrum is included in the quasar catalogue, 
while any objects originally targeted as a quasar whose spec- 
trum shows otherwise is excluded. However, quasars miscat- 
egorized likely have no spectroscopic data, and therefore will 
be missing from the QSO catalogue. 

The SDSS catalogue contains a large number of 
'serendipitous' quasars, selected often for an unusual color 
or FIRST match (|Richards et al. 2002|) . As a result, SDSS 
detection is only complete for objects brighter than 19.1 in i 
band, but includes fainter objects. The sample used in this 
paper includes all SDSS objects for which masses could be 
determined, including serendipitous objects. However, be- 
cause the SEB is determined by the most luminous and 
therefore brightest quasars at each redshift, the serendipi- 
tous sample has negligible effect upon the determination of 
the SEB. 

As part of the original target selection study, Richards 
et al. (2002) compared SDSS photometric quasar selection 
to known quasar catalogues, finding that 92.7% of 2096 
known quasars, including 94.5% of the 1540 in their 'bright' 
subsample, were correctly targeted as quasars by SDSS. 
Simulations suggest that the true quasar completeness is 
closer to 90%. There is a gap in the SDSS catalogue around 
2.5 < z < 3, where the quasar locus crosses the stellar locus. 
This gap does not affect our analysis of the SEB, which is 
restricted to z < 2. 

The SEB occurs among the brightest quasars in each 
redshift bin, so saturation in SDSS might present a problem. 
Since 3C273, at a redshift of just 0.16, is a 13th- magnitude 
quasar (Schmidt 1963), we must ensure that this luminos- 
ity cutoff is not merely an artifact of saturation. Jester et 
al. (2005) compares SDSS with the Bright Quasar Survey 
(BQS K|Schmidt fe Green 1983() derived from the Palomar- 
Green survey. Of the 51 objects in BQS that lie within the 
DR3 footprint, 29 are in the SDSS DR3 catalogue, 3 are 
above i — 15.0 and therefore excluded, and the remaining 
19 have been identified as QSO candidates in DR3 photom- 
etry but were still in the spectroscopic queue at the time 
of release. These 19 have all been targeted for inclusion in 
the DR7 catalogue, and some are included in the DR5 cat- 
alogue used here. Jester et al. (2005) also show that, while 
individual quasars are variable over 20-year time-scales, sta- 
tistically BQS and SDSS show no systematic photometric 
bias with respect to each other. 

Could the SEB occur because quasars at higher 
luminosity have different properties which cause 
them to be missed by the BQS? In related work, 
Amarie & Steinhardt (2009) compared the SDSS cat- 
alogue with the Veron-Cetty/Veron (VCV) catalogue 
( |Veron-Cetty fc Veron 2006[ l, which is a compilation of 
quasars from every available source. VCV would be a poor 



choice for a full completeness study because there is no 
guarantee of either high observational quality or statistical 
completeness in any region of the sky. However, as VCV 
includes quasars selected by all current methods, any 
population that might cross the SEB will show up. Amarie 
& Steinhardt (2009) conclude that the VCV and SDSS loci 
are well-matched, outside of the z — 2.5 — 3.0 band affected 
by the decrease in selection efficiency. Also, while VCV 
contains tens of quasars with i < 16 at low redshift, no 
substantial population of SEB-crossing quasars are found. 
We can therefore conclude that the SDSS quasar sample is 
representative of quasars as a whole, and find no evidence 
that these two new luminosity bounds are introduced by 
artificial SDSS selection. 



4.3 Bolometric luminosity errors 

The SEB takes the form of a paucity of quasars at high lumi- 
nosity. So, if the luminosity of the most luminous quasars in 
each redshift bin were underestimated, this could simulate 
such a boundary. The luminosity calculation has three prin- 
cipal components (Richards ct al. 2006b, S hen et al. 2008). 
First, the apparent magnitude in five colour bands is cal- 
culated from the SDSS photometry as part of the standard 
SDSS pipeline. Then, a K correction is used to account to 
convert the apparent magnitudes at different redshift to a 
magnitude in i band if the quasar were at z = 2. Finally, 
Richards et al. (2006) calibrate a bolometric correction to 
Mi(z = 2). Amarie et al. (2009) consider the components 
in determining the bolometric luminosity for a given quasar 
and whether such a bias might be introduced. They confirm 
that SDSS bolometric luminosities have statistical uncer- 
tainties of only a few percent for bright quasars. Systematic 
uncertainties in bolometric corrections may be far larger. 
However, explaining the SEB would require two peculiar 
properties of such a correction: (1) a systematic bias to- 
wards underestimating luminosities at high mass; and (2) a 
systematic bias with redshift that exactly cancels the mass 
bias such that the lowest-mass quasars at every redshift can 
reach their Eddington luminosity, while no quasars are ever 
super-Eddington. This combination is sufficiently improba- 
ble that we should consider systematic errors in bolometric 
luminosity estimation a very unlikely explanation for the 
SEB. 



4.4 Spectroscopic errors 

Virial mass estimates take the form given by eq. ([2|for the 
full width half maximum (FWHM) and continuum luminos- 
ity (L\) of different spectral lines. Uncertainties in these 
measurements become uncertainties in the corresponding 
mass estimate. If errors in the FWHM and L\ axe the 
source of the SEB, they must be systematic errors that lower 
L/L<Edd by leading to artificially high mass estimates, i.e., 
overestimates of the FWHM, L\, or both. 

Continuum luminosities would seem easy to mea- 
sure for SDSS spectra, with pixels of width 70 km/s 
(|York et al. 2000")) . In practice, continuum fitting is straight- 
forward near the Hj3 and CIV lines, but when using Mgll, 
if the red end of the emission line is near the 9200 A upper 
limit of the SDSS spectrograph (i.e., for z ~ 2), it may be 
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difficult to discern the extent of the iron 'bump' / Balmer 
continuum (cf. Wills et al. 1985), which must be subtracted 
as part of the continuum fit. However, even a factor of two 
error in the continuum fit would lead to just a 0.14 dex error 
in virial mass estimation. 

Fitting for the FWHM of spectral lines can be more 
difficult. There are often as few as 10-15 pixels in an SDSS 
Mgll line profile. The FWHM in the SDSS pipeline relies 
on a fit to the second moment of the line l|Shen et al. 2008). 
Some of the pixels may have poor signal-to-noise ratios or 
may coincide with an absorption line on the blue side of the 
line and should be discarded. 

To assess these affects, we performed an independent 
fit for a subsample of 3167 SDSS quasars. We fit the broad 
line shape as the sum of two Gaussian components: a broad 
component and a narrow line component (|Hao et al. 2005). 
This technique is subject to a different set of biases than that 
of Shen et al. (2008), most importantly the possibility that 
spectral lines might be skewed or non-Gaussian. Comparing 
these two techniques with different biases and sources of 
error gives a measure of systematic uncertainty in emission 
line fitting. We performed these fits on a subset of the H/3 
mass sample, where the SEB is indicated most strongly. 

Our fitting routine is based upon LINEBACK- 
FIT, part of the IDLSPEC2D utilities package 
( |Burles fc Schlegel 2006[ ). For the H/3 line, we consider 
the spectrum between rest wavelengths 4400-5200 A and fit 
a fifth-order polynomial continuum in addition to Gaussian 
[OIII] line profiles at 4959 A and 5007A as well as the 
double-Gaussian H/3 profile. Virial mass estimation requires 
the best-fitting continuum at 5100 A, so mass estimates 
are insensitive to [OIII] fitting. For the H/3 line, where the 
continuum can be well fit on both sides of the spectral line, 
these line fits were very robust with regard to individual 
pixels. This rudimentary continuum fit was incapable of 
fitting the Fell lines near Mgll, and therefore would have 
produced poor mass estimates for our Mgll sample. 

Simulating an absorption line by setting the flux in one 
pixel to zero somewhere in the spectral line resulted in virial 
mass changes of A Mbh < 0.1 dex for even the narrowest 
H/3 lines in our sample (where one erroneous pixel repre- 
sents the greatest fraction of the total line). This approach 
also appears to be only weakly sensitive to the separation 
between the broad and narrow components, and to the con- 
tinuum fit in the region surrounding the line peak. 

In Figure [121 we show the correlation between our mass 
estimates (JVfgg) and those of Shen et al. (2008; A^s nen )- 
There is a slight offset between these two mass estimates. 
However, these two virial mass estimates based upon dif- 
ferent line fitting techniques are still well-correlated, with a 
standard deviation of 0.30 dex in Mg^^/AI^ for H/3. We 
conclude that the methods used to fit the line parameters 
entered into virial mass estimates do not substantially in- 
crease virial mass estimate uncertainty beyond the claimed 
~ 0.4 dex. Moreover, as shown in Figure[9l the H/3 and Mgll 
views of the SEB at 0.6 < z < 0.8 are in strong agreement. 

The 1O 9 M0 quasars at 0.2 < z < 0.4 have a maximum 
luminosity of L/Le ~ 0.1 (Figure [3]). In order for the SEB 
to be spurious, these measurement errors would need to be 
at least 1 dex at high mass, as well as have a systematic bias 
towards overestimating the mass. The black dashed line in 
Figure [12] is drawn with the relation that would be required 
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Figure 12. A comparison of virial mass estimates from H/3 lines 
using SDSS line fits (Shen et al. 2008) and Gaussian line fits from 
this paper. The standard deviations is 0.30 dex. The red dashed 
line is drawn where the two mass estimates are in agreement. The 
black dashed line is drawn at the systematic mass offset required 
to move the SEB back to LEdd at 0.2 < z < 0.4. The number 
density has been normalized to the peak number density. 

to move the SEB back to the Eddington luminosity (Figure 
[3} • This relation is clearly larger than the differences induced 
by these two different fitting methods. We therefore find no 
evidence that measurement uncertainties in the virial mass 
estimates are introducing non-physical limits on the SDSS 
quasar locus. 

4.5 Virial mass estimates 

If the virial mass estimates are incorrect then the quasar 
locus in the M — L plane would be shifted. However, in 
order for errors in virial mass estimates to explain the SEB 
requires that the masses are incorrect in such a way as to 
produce Figure [3] by skewing the locus such that quasars 
at high mass reach either their Eddington luminosity or the 
SDSS bright-object cutoff. 

In particular, the most luminous quasars at 0.2 < z < 
0.4 have L ~ 10 46 1 erg/s (Figure [3}- In order to correspond 
to the Eddington luminosity, no quasars in this redshift 
range could have a black hole mass larger than 10 8 Mq. 
This would imply that a majority of quasars have vastly 
overestimated masses, including many central black holes 
with masses overestimated by over a factor of 10 and one by 
at least a factor of 100. Similarly (Figure [8}, at a redshift of 
1.0 < z < 1.2, the most luminous quasars have L ~ 10 47 ' 2 
erg/s, implying Mbh < 10 91 Mq. This would correspond to 
typical overestimates of at least 0.5 dex and for some objects 
0.8 dex. 

A comparison of virial mass estimates to reverberation 
mapping-based masses shows a typical uncertainty closer to 
0.4 dex (Vcstcrgaard & Peterson 2006). The virial mass es- 
timates must then have additional imprecision. Vestergaard 
& Peterson (2006) suggest these are primarily statistical in 
nature. If this is true, then if many masses in Figure 3 are 
overestimated by a factor of ten, we should also see masses 
underestimated by a factor of ten and therefore objects that 
appear to be at L = 10Ls. Since we do not, using virial mass 
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Figure 13. Quasar number density using H/3 masses at 0.6 < z < 
0.8 (red) and Mgll masses at 1.0 < z < 1.2 (blue) as a function 
of Eddington ratio for quasars within one magnitude of the SDSS 
detection limit. The dashed line is an exponential decay with an 
c- folding of 0.15 dex in L/Le- 

mis-estimation as an explanation for these new boundaries 
would require a systematic component as well. We attempt 
to evaluate both statistical and systematic uncertainties in 
virial mass estimates below. 

4-5.1 Virial mass estimates: statistics 

We can make an independent estimate of the statistical error 
in Mbh by considering the low- mass end of the M — L plane 
at each redshift, where some quasars reach their Eddington 
luminosity. Statistical uncertainties should result in some 
quasars lying at a lower mass than is physically allowed and 
thus appearing to be above LecLcL- We can estimate the true 
statistical uncertainty in Mbh from the width of this falloff 
in quasar number density. 

The quasar populations within one i-band magnitude of 
the SDSS detection limit (Figure [T3)) do not reach LEdd, but 
decline sharply at high Eddington ratio, with an e-folding 
number density rate of ~ 0.15 dex in L/LEdd- These quasars 
approach LEdd more closely than at higher luminosity, where 
quasars are bounded by the SEB. This ~ 0.15 dex falloff is 
a combination of statistical uncertainty (which should be 
Gaussian), possible systematic effects, and possibly a real, 
underlying decline in the quasar L j 'LEdd distribution associ- 
ated with the Eddington limit. Therefore, without knowing 
the details of underlying physical cause of the decline, there 
is a ~ 0.15 dex upper limit for the statistical uncertainty in 
the Shen et al. (2008) virial mass estimates. If the underly- 
ing physical cause is a gradual rather than a sharp decline 
in number density, the maximum possible statistical uncer- 
tainty in virial mass estimates will be reduced. 

0.15 dex is notably smaller than the uncertainty esti- 
mated by Vestergaard & Peterson (2006) in a comparison of 
reverberation and H/3-based virial masses for the same set 
of objects. This suggests that the disagreements between 



the reverberation mapping and virial mass estimates have a 
substantial systematic component. A closer examination of 
the Vestergaard & Peterson (2006) mass estimate compari- 
son suggests that reverberation mapping-based estimates at 
high mass might be larger than virial mass estimates. This 
would result in a stronger SEB (i.e., one with a lower slope) 
than reported in Table [4] 

The statistical uncertainty in Mbh can also be esti- 
mated by comparing different virial mass scaling laws at 
redshifts where two lines can be used. At 0.6 < z < 0.8, 
there are 3505 quasars in common between the H/3 and Mgll 
samples. Since the SEB occurs at the high-luminosity end 
of the quasar sample, if is due to statistical uncertainties in 
mass estimation there should be a higher dispersion between 
different mass estimates for more luminous quasars than for 
the population as a whole. 

At 0.6 < z < 0.8, the mass dispersion for all quasars 
common to the H/3 and Mgll samples is 0.20 dex, while the 
mass dispersion for the 10% most luminous quasars is 0.16 
dex. These dispersions are substantially smaller than the ~ 
1.0 dex required to remove the SEB. The smaller dispersion 
for luminous quasars likely occurs because brighter objects 
are typically accompanied by higher signal-to-noise spectra. 
Not only are individual H/3 and Mgll-based mass estimates 
well-correlated for bright quasars, but so is the resulting 
SEB. At 0.6 < z < 0.8, the best-fitting SEB shown in Table 
rjhas slope a = 0.60±0.06 for H/3 masses and a = 0.61±0.10 
for Mgll masses. 

Statistical uncertainty could produce high-mass outliers 
that might contribute to high-mass deviations from the best- 
fitting line, but a correction of the SEB to LEdd would re- 
quire an overestimate of at least 1.0 dex for many quasars. 
These estimates for the statistical uncertainty in Mbh range 
from just 0.15-0.4 dex, too small to artificially produce the 
SEB. Further, statistical uncertainty would broaden the lo- 
cus, not change the slope, and the SEB has a slope at least 
3.7<7 away from the a = 1 of LEdd in every redshift bin. We 
conclude that while statistical uncertainty is a factor in the 
final shape of the quasar locus at the 0.15 dex level, we find 
no evidence from any of our tests that statistical uncertainty 
could be responsible for the ~ 1 dex SEB shown in Figure 

El 

4-5.2 Virial mass estimates: systematics 

The falloff in Figure [T3] starts at log LjLsdd ~ —0.5 at low 
mass and log L/LEdd ~ —1.0 at high mass for H/3. If some 
quasars reach their Eddington luminosities, this early falloff 
suggests that masses could be systematically overestimated. 
This might mean masses are overestimated by as much as 
0.5 — 1.0 dex. A systematic overestimate in virial mass esti- 
mates would likely also require a systematic overestimate in 
the reverberation mapping-based estimates against which 
they are calibrated, although a sufficiently small overesti- 
mate might lie within the 0.4 dex uncertainties between the 
two methods. It appears that systematic errors of of ~ 1 dex 
in virial mass estimation cannot be entirely ruled out with 
this analysis. 

However, the SEB slope requires not just a shift in Mbh 
but a systematic, mass-dependent change in the mass es- 
timate. Low-mass quasars cannot have mass overestimates 
large enough for their luminosities to exceed LEdd, so that a 
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correction of the SEB might require a 1.5 dex overestimate 
at the high-mass end with only a 0.5 dex overestimate at 
low mass. Since luminosities at 0.2 < z < 0.4 (Figure O run 
over a range of only 1.3 dex, the actual underlying black hole 
masses must run over a range no larger than 1.3 dex them- 
selves, or else the SEB cannot really have slope a = 1. If 
true, this would replace the SEB with a different surprising 
feature of the quasar M — L distribution. 

For Mgll masses in particular, Onken & Kollmeier 
(2008) examine SMBH for which both H/3 and Mgll masses 
are available and, assuming the H/3 mass is the better indica- 
tor, find that the Mgll-based Mbh may be overestimated at 
high Eddington ratio and underestimated at low Eddington 
ratio. Correcting for this effect will drive low-mass objects 
at high Eddington ratios to even lower masses and there- 
fore closer to L B dd, but high mass objects at low Eddington 
ratios to higher masses and lower Eddington ratios. There- 
fore, such a correction could produce a stronger, more sub- 
Eddington boundary. 

Risaliti, Young, & Elvis (2009) also compare H/3 and 
Mgll masses, finding that 

\og[M BH (Hp)] = 1.8 x \o g [M BH {MgII)] - 6.8. (3) 

This correction acts in the same direction as the correc- 
tion proposed by Onken & Kollmeieier and will change the 
SEB slope a in each Mgll mass bin from L — aM to 
L — a' (1.8M'), for a reduction in a by a factor of 1.8 (Figure 
[9]). This correction produces a 3<r disagreement between the 
H^ and Mgll SEB slopes at 0.6 < z < 0.8, although since 
these two samples comprise a very similar set of objects, it 
is almost certainly incorrect to treat their SEB uncertain- 
ties as statistical in nature. This result seems to require that 
the brightest quasars at each mass need smaller Mgll mass 
corrections than fainter ones. Perhaps high signal-to-noise 
spectra yield more reliable Mgll masses. We conclude that 
specific, known systematic errors in virial mass estimation 
might take the SEB slopes closest to Eddington and move 
them further sub-Eddington. Any new correction that might 
bring SEB slopes to Lsdd would replace the SEB with a dif- 
ferent surprising feature of the quasar distribution in the 
M — L plane. 



4.6 Underlying mass distribution 

One may wonder whether we have simply re-established the 
well-known decline of the luminosity function of quasars at 
high luminosity (Schmidt 1968. Richards et al. 2006a), and 
whether the 'missing' objects simply reflect the underlying 
mass distribution. In Figure [14] we compare the mass dis- 
tribution of black holes in two well-populated luminosity 
bins at the same redshift (0.2 < z < 0.4) separated by 0.2 
dex. If the SEB were the expression of mass turnoff, objects 
in the bin 0.2 dex higher in luminosity would also be 0.2 
dex more massive, just that there would be fewer at higher 
mass and luminosity. In fact, the higher-luminosity objects 
are on average just 0.09 dex more massive. Further, a KS 
test comparing the mass distributions with a shift of 0.09 
dex results in a probability of 47.2% that they are identical, 
while a shift of 0.2 dex results in a probability of only 6.0% 
that the two are drawn from the same distribution. Not only 
does the SEB have a slope a < 1, but for the entire quasar 



,11 




_n 










1 1 1 1 1 1 


: J 

, ft 












_r 

u 


_TL 


7 


8 












9 10 



Log(M/M e ) 

Figure 14. The quasar mass distribution for two different lumi- 
nosity bins in the 0.2 < z < 0.4 sample shown in Figure [3] We 
divide the sample by luminosity into bins with 10 45 ' 3 < L/Lq < 
10 45 ' 5 (601 objects, red) and 10 45 5 < L/Lq < 10 45 ' 7 (203 ob- 
jects, green). The higher-luminosity objects are on average 0.20 
dex higher in luminosity but only 0.09 dex higher in mass. 



population, an increase in mass corresponds to a sub-linear 
increase in luminosity. 

In particular, the SEB is not simply the expres- 
sion of mass turnoff, as the luminosity increase is only 
slightly correlated with mass increase at the bright end. 
This does raise the possibility, discussed in Paper II 
(Steinhardt & Elvis 2009), that we can probe quasar turnoff 
using SDSS data. 

The SEB is an expression of the relative scarcity of high- 
mass, high-Eddington ratio objects compared to high-mass, 
low-Eddington ratio and low-mass, high-Eddington ratio ob- 
jects at each redshift. While the mass estimate does depend 
upon luminosity (as well as the FWHM of a broad emission 
line), an increase of 1 dex in luminosity leads directly to an 
increase of just 0.5 dex in mass (Eq. [2}, not the ~ 1.7 dex 
increase required for a typical SEB slope. Since the SEB is 
determined using objects only at a common redshift, there 
are no changes in the comoving volume to account for in 
evaluating this scarcity. 

Because both masses and (bolometric) luminosities are 
derived quantities, the SEB could be alternatively inter- 
preted as either (1) a physical clue about the components of 
these derived quantities (e.g., at high masses, an increase in 
luminosity might in every case be accompanied by a much 
larger than expected increase in the FWHM of broad emis- 
sion lines), (2) a failure in deriving these quantities (e.g., 
virial mass estimation breaks down on the high-mass end 
at each redshift in exactly such a manner as to produce 
this skew but is correct at higher redshift when identical 
masses lie on the low-mass end), or (3) a true lack of high- 
mass, high-luminosity quasars. We have argued that (3) 
should be the favored interpretation because of the fine- 
tuning required if (1) or (2) is to apply only the highest- 
mass objects at each redshift. However, all three interpre- 
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Table 6. Possible explanations for the sub-Eddington boundary 
considered in § |4] 



Explanation 


Plausible? 


Statistical Insignificance 


No 


SDSS Selection 


No 


Luminosity Errors 


No 


Line Measurement Errors 


No 


Virial Masses: Statistics 


No (maybe at high M) 


Virial Masses: Systematics 


Unlikely 


Quasar Turnoff 


No 



tations would require that something physical is different 
about high-mass, high-luminosity quasars in each cosmolog- 
ical epoch and therefore any interpretation presents a theo- 
retical puzzle. 

We have considered six potential explanations for these 
new limits without introducing new quasar physics, as sum- 
marized in Table [6] A large, unknown systematic error in 
virial MbH measurements is the most plausible explanation 
by comparison for how the SEB might be produced without 
additional underlying physics, but requires that virial mass 
estimates are ~ 1 dex incorrect, reverberation mapping is 
~ 1 dex incorrect, and that quasars at low redshift lie within 
a spread of just 1.3 decades of central black hole mass. It 
would appear likely that the explanation for these new limits 
lies in new theory rather than observational errors. 



5 DISCUSSION 

Quasar catalogues such as the Sloan Digital Sky Survey are 
now large enough to yield useful statistical conclusions after 
being subdivided. In this work, we have explored just a few 
of the many ways in which the catalogue might be divided. It 
is immediately apparent that when the quasar distribution 
is considered at any given redshift, the most massive central 
black holes do not accrete at their Eddington luminosity, 
but rather all fall well short of Lb. While there are many 
measurements that contribute to these quasar distributions, 
each with its own set of assumptions, a close examination re- 
veals no evidence that any these are responsible for massive 
quasars remaining sub-Eddington. Rather, it appears that 
this 'sub-Eddington boundary' is a new physical limit for 
quasars. Thus, quasar accretion must be more complicated 
than had been previously thought. 

The sub-Eddington boundary helps to recast the prob- 
lem in a two-dimensional form. This form emphasizes the 
luminosity as being a function of the mass and accretion 
rate, fundamental properties of the quasar, in a non-trivial 
way. Quasar luminosity functions and quasar mass functions 
are both one-dimensional projections of the mass-luminosity 
plane. The full two-dimensional quasar distribution contains 
complexities difficult to discern from either of these projec- 
tions. 

The sub-Eddington boundary presents a theoretical 
puzzle. Not only must an improved theory explain the sub- 
Eddington boundary, but it must also explain its detailed 
shape and explain the evolution of that shape with red- 
shift. While the slope of the sub-Eddington boundary may 
not vary greatly with redshift, the mass scale at which the 



boundary becomes relevant is larger at higher redshift. The 
results in this paper can be used to develop a series of tests 
which should be capable of discriminating between models. 

Given the complications introduced by the sub- 
Eddington boundary, the case could be made that we re- 
ally understand surprisingly little about the supermassive 
black holes that appear to be at the center of nearly every 
galaxy. Not only is their seeding mechanism unknown, but as 
shown in this paper, the growth mechanism is quite poorly 
understood, contrary to expectations. The Eddington limit 
is relevant at low black hole masses, but is only part of the 
story. The sub-Eddington boundary developed in this work 
is the latest addition to a growing collection of puzzles re- 
garding every phase of the evolution of galactic nuclei and 
surrounding regions. 
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